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THE DISTRIBUTION OF GENE FREQUENCIES UNDER 
IRREVERSIBLE MUTATION 


By SEWALL WRIGHT 
DEPARTMENT OF ZOOLOGY, THE UNIVERSITY OF CHICAGO 


Communicated June 13, 1938 


Under reversible mutation, the frequency (q) of a gene, subject to sys- 
tematic evolutionary pressure (Ag) and to the accidents of sampling in a 
limited population (N diploid individuals), varies according to a certain 
distribution (¢(q)) discussed in previous papers.'?* If mutation is irre- 
versible, the distribution curve for such genes should attain constancy of 
form, but all class frequencies should fall off at a uniform rate (K) as genes 
drift irreversibly into fixation. The purpose of the present paper is to 
broaden somewhat the treatment in this latter case. 

As previously shown? the rate of fixation is approximately half 
the frequency, f(l1 — 1/2N), in the subterminal class. Thus with 


= 1 


(1 — 1/2N) 


K = 1/2f(1 — 1/2N) = © aN approximately. (1) 


1 
The changes in the mean (g = - g¢(q)dq) and the variance (6, = 
0 


1 
[ (q — G)*¢(g)dq) of gene frequencies, due to fixation in one generation, 
/ 0 


may be expressed as follows in terms of the systematic evolutionary pres- 
sure, Ag, and the variation due to sampling, ong 


Aq — + = (1 — + (8) 


The latter can be reduced to following, ignoring a negligible term in 
(Ag)? 


= 


254 GENETICS: S. WRIGHT Proc. N. A. S. 


2 f (q — + o4,¢(q)dq = K[(1 — g)* — (4) 


If the conditions are such that under equilibrium, with reverse mutation 
at an indefinitely low rate, there is no important accumulation of genes in 
the class g = 1, complete irreversibility of mutation should make no 
appreciable difference in the form of distribution. The demonstration of 
the formula for this case (K = 0) can be put in a very simple form. 


Let Agqe(g)dq = x(q). (5) 
Equations (2) and (4) can be written as follows, putting K = 0, 
x(1) — x(0) = 0. (6) 
f xa txt) + 0 Dx) = 5 f 0. 
Equation (7) is obviously satisfied by the following 
x(@) — fax) + (1 =0. 8) 


This means little, until it is shown that (8) also satisfies (6). 

As there can be no sampling variance in homallelic populations, o4, = 0 
ifg = Oorg = 1. Equation (8) reduces to (6) if g = 0 org = 1 and both 
¢(0) and ¢(1) are finite. Equation (8), therefore, satisfies the condition of 
constancy of the mean as well as that of constancy of the variance. 


[x(q) — x(1)] can be evaluated as follows, using (5) and (9), 


dx(q) 2 Aqe(q)dq 
(12) 
Taq 


This is the desired expression for the distribution in terms of Ag and 


2 
Taq . 


Putting o4, = , its value in a population of N diploid indi- 


me 
2N 


_ 2Ix(@) = (9) 
= 
o Aq 
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viduals, it is the same as the formula given previously* (except for the 
value of the coefficient C). 

If K is not zero, this method does not appear to lead to a usable general 
expression which satisfies (1), (2) and (4). An expression for the rate of 
decay can, however, be obtained from (2). Let v be the rate of mutation, 
and let W be the mean selective value of genotypes. As shown previously* 


dlog W 


Aq = — g) + — (13) 


— — ge(g)d log W. (14) 
2(1 0 


If there is no selection (W = 1) 
K =v. (15) 


The distribution (16) for this case was obtained in a previous paper? by 
substitution of Ag = v(1 — gq) in what is equation (18) of the present 
paper. It can easily be verified that (16) also satisfies (1), (2) and (4). 


= (16) 


The most important cases are probably those in which the size of popu- 
lation is so small that recurrent mutation has no effect on the form of 
distribution (v much less than 1/4N). The selection pressure for any 
degree of dominance can be written sufficiently accurately 


Aq = (s + ¢g)q(1 — q). (17) 


The condition that the frequency of any class of gene frequencies, q,, 
be reconstructed after each generation, except for a uniform decay at 
rate K can be represented as follows,” using p = 1 — gq for brevity. 


(1 — K)¢(q.) = A (q + — 


r(2N) 
~ pagel (2Na.) (18) 


where A 


If v is so small that practically all genes are fixed in the class f(0), K 
may be ignored in determining the form of the distribution for unfixed 
genes. Substituting the value of 4g from (17) 


o(q.) =A f + p(s + tg) [1—g(s + tg) (19) 


The following approximations may be used 


[1 + pls + tg) = [1 — Ng.p%(s + (20) 
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[1 — g(s + tg) = TL — Npgr(s + (21) 
The product of expressions (20) and (21) is approximately 


Since the random deviations of g have the variance, = — 


is of the order 1/2N. The term N(s + tq)*(q, — q)? is thus negligibly 
small compared with N(s + tq)*p,q, which itself is as small a term as it is 
necessary to consider. The former may thus be ignored. The constant 
and variable gene frequencies are separable in the exponential term in (22) 
except in the term e~*"~®*, — The exponent in this case is smaller than 
—tand the term can be written [1 — N¢(q, — q)?] with sufficient accuracy. 
Equation (19) can now be written as follows: 


0 
Nt(q. — 9)*]¢(g)dg. (28) 


(Co + Cig + Caq?...). (24) 


2Nsq + Nig? 
q (1-4) 


This entirely eliminates the exponential terms, leaving (23) in a form 
which can be solved, using the following sufficiently accurate formula in 


Let eq) = ° 


which terms of the order - are ignored. 


I 


The resulting coefficients of the powers of g, on the right side may be 
equated to those on the left side leading to the following general expression 
(in which = = C-3 = 0). 


m(m + 1) ce. 


Cn—2 (2Nst)—Cy—3 Nt. (26) 


The higher coefficients can all be expressed in terms of Cy and C, and 
substituted in (24). By letting C, = 2NsC, + D the terms for which C) 
is the coefficient can be condensed into exponential form. It will be con- 
venient to substitute C for Cy. 


2N. Nac 2N. Ntq? 
t — Nb 
x(x 1) 25 
Cm 4N 2 
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Nsa + Ntq? on 
og) = + Davia) | (27) 
where 
4 
1 (2Nsq) 2(2Nsq)? , 2(2Nsq)* , 3(2Nsq)4 
4 4 2(2Nsq) 69(2Nsq)? she ) 
5 (5 [7 
+ (2...) + 
IfD =0 
Cet + 2Ntq? 
29 
o@) = (29) 


This is the case of equilibrium under reversible mutation or irreversible 
mutation opposed by sufficiently strong selection as can be seen by sub- 
stituting Ag = (s + ég)g(1 — Q), 
except in the coefficient. 

The case of irreversible mutation with fixation occurring at a low rate, 
can be found from (1) and (2), assuming that nearly all genes are in one 
of the homallelic classes. It should be noted that the formula for (gq) 
only applies where K is of lower order than 1/2N, 2Ns* + tin (26). Muta- 


tions to the class g = 1/2N contribute the amount 2Nzf(0) = fa /2N) 


and these contribute to the change of mean by the amount f(1/2N)/4N. 
The mean, however, must be so low that the term K(1 — g) may be 
written K sufficiently accurately. The following relations are all approxi- 


in (12). They agree 


mate: 
1 1/2N (1 — 1/2N) 
a 1 D 
f(1/2N) = chi 2s | +2. (31) 
f(l 1/2N) = 
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1 4Ns + 2Nt 1 Ds Dt 
q¢(q)dq | aN an | + 5 (33) 


Substituting (31), (32), (33) in (30) and substituting y(1) for y(1 — 1/2N) 


2Ns? 
(leading terms in difference, t/3, 3 ) 
(4) 
With practically all genes in the class g = 0, f(0) = 1 = —“—~= -—_. 
ith practically all gene q f(0) a: “ie 


Thus C = 4Nv approximately. 


4Nve*N + er + (35) 
q(1 — q) ¥(1) 


For sufficiently small values of Ns and Nt we may take ¥(g) = 1 + 
1/3Ntg? and represent the exponentials by the first two terms of their 
expansions. The following shows how the hyperbolic distribution 4Nv/q 
is modified by weak selection (s positive for favorable mutation, negative 
for unfavorable mutation). 


= 


o(q) = [1 + 2Nsq + 1)]. (36) 


The rate of fixation (K) of genes can be found from the left member of 
(30) (in which inaccuracies in the evaluation of D have less effect than in 
the right member). 

¥(1) 

This reduces to v(1 + 2Ns + 2/3Nt) for such small values of Ns and 
Nt as implied in (36). It is to be noted that irreversible mutation should 
ultimately lead to fixation of the mutant even when opposed by selection 
(s negative) but the rate is exceedingly slow unless Ns and Nt are small. 

In the special case of genic selection (¢ = 0) 

2Nsq 2Nsa 


— 
¥(q) = (38) 
4Nv 
(39) 

An essentially similar derivation of this formula has been given previously 

4Nvs 
by the author’ and a different one by Fisher.‘ In this case K = oo 


approaching v(1 + 2Ns) as 4Ns decreases. 
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The question of the chance of fixation of an individual mutation must 
be distinguished from the rate of fixation (K) under recurrent mutation. 
The chance of fixation is given by the ratio f(1 — 1/2N)/f(1/2N) = K/2Nv. 
In the case of no dominance, this gives 2s/(1 — e~*’*) or approximately 
2s for favorable mutations occurring in a large population, in agreement 
with Fisher.‘ For indifferent factors it is 1/2N Unfavorable mutations 
have a chance of fixation 2s/(e*’’ — 1) but this is small unless 4Ns is 
small. 

The results presented here bear on the possibility of a course of evolu- 
tionary change determined by mutation pressure, a process which at first 
sight seems the most obvious implication of modern genetics. The pos- 
sibility does indeed exist but requires either an almost complete indifference 
of the mutation with respect to adaptive value or else a very small effective 
size of population over a long period of time. The most important case in 
which mutation pressure seems likely to be a major factor is that of extreme 
degeneration or elimination of organs that have ceased to be useful.®® 
The degeneration of the eyes and loss of pigment of cave forms is an 
example of a case in which the conditions make it especially probable that 
mutation pressure is a real factor. In all of these cases, however, the 
likelihood that various direct and indirect effects of selection may also 
play a réle should not be ignored.® 

It should be noted that while the average rate of fixation of irreversible 
mutations is low, the large element of chance with respect to which muta- 
tions become fixed in each particular case makes this a greater factor in 
the diversification of small isolated populations than is at first apparent. 
Indeed there may be much diversification of gene frequencies among such 
populations under conditions in which there is no appreciable systematic 
tendency toward fixation of the sort investigated here. 


1 Wright, S., Amer. Nat., 63, 556-561 (1929). 

2 Wright, S., Genetics, 16, 97-159 (1931). 

3 Wright, S., Proc. Nat. Acad. Sci., 23, 307-320 (1937). 

4 Fisher, R. A., Proc. Roy. Soc. Edinburgh, 50, 205-220 (1930) 
5 Wright, S., Amer. Nat., 63, 274-279 (1929). 

6 Haldane, J. B. S., Jbid., 67, 5-19 (1933). 
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AUXIN IN ISOLATED ROOTS GROWING IN VITRO 
By J. VAN OVERBEEK AND JAMES BONNER 


Wiiiiam G. KERCKHOFF LABORATORIES OF THE BIOLOGICAL SCIENCES, CALIFORNIA 
INSTITUTE OF TECHNOLOGY 


Communicated June 15, 1937 


Up to the year 1933 opinion was divided as to whether or not auxin is 
present in the root. In this year Boysen-Jensen! conclusively demon- 
strated that auxin is present in the normal root. With respect to the iso- 
lated root, however, Fiedler* has reported that auxin, although initially 
present, disappears completely and in general within 24 hours, when the 
root is cultivated im vitro. Nagao* on the other hand, has shown that 
auxin may be recovered from isolated roots after 6 days’ cultivation in 
vitro. The object of the present work has been primarily to establish, 
with the aid of improved techniques, whether or not auxin is actually 
present in isolated roots cultivated in vitro. A second but related problem 
of whether or not auxin is produced by such roots will be dealt with in a 
later communication although experiments of a preliminary nature indi- 
cate that it is not produced but is merely carried along in the tip of the 
isolated root. 

Methods.—Pea roots were cultivated in vitro under the conditions which 
have been described elsewhere.‘ Tips ca. 4 mm. long were removed from 
the roots of germinating pea seeds® and cultured under strictly sterile con- 
ditions in 10 cm. Petri dishes in a basic nutrient solution containing sucrose, 
inorganic salts, vitamin B, (0.1 mg. per liter) and varying amounts of 
certain amino acids. After the roots had grown for one week and were ca. 
70 mm. long, they were subcultured by removal of a 10 mm. tip to fresh 
medium. If the roots were to be further subcultured this procedure was 
repeated at weekly intervals. Since lateral roots are formed only after 
approximately 10 days’ culture im vitro under these conditions, these ex- 
periments are uncomplicated by the presence of more than one growing 
point per root. 

The technique of the auxin extraction and determination has been de- 
scribed in detail elsewhere.° The material to be extracted was placed in 
highly purified ether, in the cold room (ca. 2°C.) for 15 to 30 hours. The 
material was neither ground nor acidified’ since either of these procedures 
may result in a partial destruction of auxin. The ether extract was then 
evaporated to dryness and taken up in a known amount of hot agar (1.5%). 
This agar was in turn analyzed for auxin by means of the Avena test. 
For each determination 24 test plants were used. In order that the 
auxin content of the roots might be expressed in ‘‘indole-acetic acid equiva- 
lents,’’ comparable controls with this substance were also run. In general 


i 
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two concentrations, one of 20 gamma and one of 7 gamma per liter were 
used for these controls. Each control was also run with 24 test plants. 
A blank with plain agar blocks was also included in each test. In an 
earlier paper it has been pointed out that the Avena curvature-concentra- 
tion curve often does not pass through the origin, but that it intersects the 
abscissa. This means that the curvature is not directly proportional to the 
concentration. There is a certain “threshold concentration’”’ before a 
curvature is obtained. This threshold concentration, which is in general 
about 5 gamma per liter, was determined for most of the tests and was ap- 
plied to the subsequent calculations. It was, however, found that if this 
threshold concentration is taken into account, the relationship between 
auxin concentration and Avena curvature can be expressed by a straight 
line,’ at least under the conditions of these experiments. 

The auxin content of the roots was calculated according to the following 
formula: 


(C x + O) x 
WwW 


auxin concentration = 


gamma indole acetic acid 
equivalents per kg. fresh 
weight 


the average curvature obtained in the Avena test (extracted 
material) 

indole acetic acid concentration giving an increase of 1° in the 
Avena test, in gamma per liter 

threshold concentration, gamma indole acetic acid per liter 

volume of agar in which the dry extract was taken up, cc. 

= fresh weight of material under investigation, gms. 


ll 


II 


Auxin Concentration in Roots.—Table 1 shows that under the conditions 
of our experiments auxin is present in isolated pea roots cultivated in vitro 
even after three weeks, and consequently, after 2 subcultures. There can 
then be no doubt but that auxin may be present in relatively large amounts, 
even in the cultured root. Fiedler,? as mentioned above, failed to find 
auxin in isolated roots (Zea, Pisum, Vicia) after more than 48 hours’ 
growth in vitro. That Fiedler failed to extract auxin from Zea roots of 
any age, is undoubtedly due to the extraction technique which he em- 
ployed. It has been shown elsewhere® that auxin cannot in general be 
obtained from Zea if acid is employed in the extraction (acid was used 
by Fiedler). It is shown in table 2 that with the present extraction tech- 
nique, auxin can be obtained from the roots of sterile, germinating Zea 
seeds. The amount present in these roots is smaller than that present in 
Pisum roots of the same length (table 2). It has been shown earlier® that 
the shoots of Zea seedlings, also, contain less auxin than do the shoots of 
corresponding Pisum seedlings. The failure of Fiedler to extract auxin 
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from isolated pea roots cannot, however, be satisfactorily explained, al- 
though this also may have been due to the different techniques employed. 

Nature of the Hormone Present.—It was of importance to ascertain 
whether the auxin of these roots was auxin-a, auxin-b or hetero-auxin 
(indole acetic acid). This may conveniently be determined by the dif- 
ferential acid-alkali destruction test of Kégl, Haagen-Smit and Erxleben. 


TABLE 1 


AvuxIn CONCENTRATIONS FouND IN ISOLATED Roots GROWING in vitro, GAMMA “IN- 
DOLE AcETIC AcID EQUIVALENTS” PER KG. FRESH WEIGHT. THE CULTURES WERE 
ANALYZED ONE WEEK AFTER THEY WERE STARTED. CULTURES OF ‘‘ALASKA”’ Roots 
ARE INDICATED BY A, THOSE OF “‘PERFECTION”’ BY P. I, II, III, Etc., ADDED TO AOR 
P INDICATE RESPECTIVELY THE FIRST CULTURE, SECOND CULTURE (FIRST SUBCULTURE), 
THIRD CULTURE, ETC. THUS THE P III CULTURE AT THE TIME IT Was ANALYZED Hap 
BEEN CONTINUED THREE WEEKS SINCE 1HE INITIAL Root Tip Was Cur FROM THE 


EMBRYO 
EXPT. NUMBER CULTURE 7/KG. REMARKS 
80315 Al 1.69 
80413 Al 2.24 no tip 
80418 All 5.88 
7.48 
80429 All 6.40 
, 80517 roots 4 mm. long from embryo, P 0 22.3 
80315 PI 2.26 
80315 PII 2.28 
80331 PII 2.45 
80418 PIII 4.12 
TABLE 2 


AUXIN CONCENTRATION AND CONTENT OF 1-cM.-LONG Roots oF GERMINATING Zea, 
AND Pisum SEEDS. GAMMA “INDOLE ACETIC ACID EQUIVALENTS” 


PLANTS C (CXh° +0) Vagzar Ww +/KG. /PER ROOT 
80603 Zea 9.4 0.5 2.60 1.64 20.6 X 10-* 


80523 Alaska 
80523 Perfection 


23.6 10.0 6.60 36.0 300.0 X 10% 


18.0 7.5 4.00 i 
9.8 15.0 re 35.3 268.0 X 10 


This test is based on the facts (a) that indole acetic acid is heat stable to 
alkali but is destroyed by acid, (5) that auxin-a is stable to acid but is de- 
stroyed by alkali and (c) that auxin-b is unstable to both acid and alkali. 
Three separate tests of this nature showed that the auxin extracted from 
pea roots was completely destroyed by alkali, unaffected by acid and must 
hence have been auxin-a. It might be mentioned in this connection that 
Heyn!! has shown by another method that the auxin of the roots of ger- 
minating Vicia seeds is also auxin-a. It might also be noted, however, 
that one experiment made with non-sterile roots of pea seedlings showed the 
presence of large amounts of indole acetic acid, presumably formed through 
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the action of microérganisms. Reserve must therefore be exercised in the 
interpretation of experiments dealing with the auxin relations of non-sterile 
roots. 

The experimental procedure was as follows: A large amount (45 gms.) 
of roots which had been one week in culture and most of which were with- 
out tips (tips had been removed for the subsequent culture) was extracted 
with ether in the usual fashion, and the extract divided into 4 equal parts. 
One part was analyzed immediately for the determination of the amount of 
auxin originally present. The other three portions were placed in flasks 
and the ether removed completely by distillation. To one sample was 
then added 3 cc. of 5% HCl; to the second sample, 3 cc. of 0.5 normal 
NaOH and to the third sample 3 cc. of distilled water. The samples were 
then refluxed on a boiling water bath for 15 minutes, cooled, neutralized 
(slightly acid) and extracted with ether. This ether was then tested in the 
manner described above. The results of one experiment are summarized 


TABLE 3 


DETERMINATION OF THE NATURE OF AUXIN IN STERILE, ISOLATED Roots GROWING 
in Vitro. Expt. NUMBER 80509 


REFLUXED 


WITH Cc (C +O) Vagar Ww 
HCl 2.6 0.6 X 44.6 0.420 
NaOH +0.9 0 0.6 1), X 44.6 0 
H,0 2.0 6.9 0.6 1/, X 44.6 0.378 
not refluxed 3:3 7.5 1.0 1/, X 44.6 0.681 


in table 3. It is of interest to note that the entire analysis was carried out 
with 30 millionths of a milligram “indole acetic acid equivalents.’ This 
corresponds, according to Kégl, Haagen-Smit and Erxleben,!? to 15 mil- 
lionths of a mg. of auxin-a. 

Distribution of Auxin in the Root.—Table 4 shows that pea roots both 
at the end of the first culture, i.e., one week in vitro, and at the end of the 
second culture, i.e., two weeks in vitro, exhibit a very pronounced gradient 
of auxin concentration from the tip toward the base. This gradient is 
much steeper than that found by Thimann’ for roots of Avena seedlings, but 
is similar to that found by Boysen-Jensen! for roots of germinating Zea 
and Vicia seeds. 

Is Auxin Produced by Isolated Pea Roots in Vitro?—This question will be 
treated in detail in a later communication. However, the present experi- 
ments permit of a comparison of the total amount of auxin per average 
initial root tip (the 4 mm. tip cut from the germinating seed) with the total 
amount of auxin per average root_after it has been two weeks in culture 
(and has hence been subcultured once). The average initial root tip (4 
mm.) contains 137 X 10-* gamma indole acetic acid equivalents. The 
total auxin content of the average root at the end of the second week is 
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TABLE 4 


DISTRIBUTION OF AUXIN-a IN STERILE, ISOLATED Roots Grown in Vitro, GAMMA 
“INDOLE AcETIC EQUIVALENTS”’ PER KG. FRESH WEIGHT 


EXPT. LENGTH 
NUMBER CULTURE PART (AVERAGE) Vagar Ww 7/KG. 
80524 ge | tip 15 mm. 11.2 25.3 0.5 1.001 12.7 

middle 15 mm. 2.1 8.9 0.5 2.003 2.22 
base 24 mm. 15.9 34.0 0.5 7.609 2.23 
80511 PII tip 15 mm. 2.6 8.7 0.6 0.856 6.10 
middle 25 mm. 19 7.5 0.4 2.500 1.20 
base 25 mm. p Ee | 6.0 0.5 8.900 0.337 
80516 PII tip 15 mm. 5.3 18.8 0.5 1.020 9.2 
middle 15 mm. 4.6 16.4 0.5 2.695 3.04 
base 20 mm. 3.2 12.0 0.5 8.490 0.71 


70 X 10~-* gamma indole acetic acid equivalents. Thus although the root at 
the end of the second week does contain auxin, it appears to contain sub- 
stantially less of this substance than does the initial tip. 

Summary.—(1) Isolated Pisum roots cultivated in vitro were found to 
contain auxin for at least three weeks after the original tip was removed 
from the germinating seeds. (2) This auxin obtained from roots under 
sterile conditions has been shown to be auxin-a. (3) A steep auxin gra- 
dient was found in these isolated roots, the highest concentration being found 
in the tip. (4) Roots after two weeks’ cultivation im vitro appear to con- 
tain less auxin than did the initial root tips. 


1 Boysen-Jensen, P., Planta, 19, 354 (1933). 

2 Fiedler, H., Zeit. Bot., 30, 385 (1936). 

3 Nagao, M., Sci. Rep. Tohoku Imp. Univ., 12, 191 (1937). 

4 Bonner, J., and Addicott, F., Bot. Gaz., 99, 144 (1937). 

5 Pea seeds of the variety ‘‘Perfection,’’ supplied by the Ferry-Morse Seed Co., San 
Francisco, were used. Seeds of the variety ‘‘Alaska” were used in a few experiments 
with substantially the same results. 

6 van Overbeek, J., Proc. Nat. Acad. Sci., 24, 42 (1938) and Bot. Gaz. (1938) (in press). 

7 Thimann, K. V., Jour. Gen. Physiol., 18, 23 (1934). 

8 yan Overbeek, J., [bid., 20, 283 (1936). 

® Compare: Thimann, K. V., and Schneider, C. L., Amer. Jour. Bot., 25, 270 (1938). 

10 Kégl, F., Haagen-Smit, A., and Erxleben, H., Zeit. Phystol. Chem., 228, 104 (1934). 

11 Heyn, A. N. J., Proc. Kon. Akad. Wetenschap., Amsterdam, 38, 1074 (1935). 

12 Kgl, F., Haagen-Smit, A., and Erxleben, H., Zeit. Physiol. Chem., 228, 90+(1934). 

13 Report of work carried out with the aid of the Works Progress Administration, 
Official Project Number, 465-03-3-342, Work Project N-9199. 
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THE MECHANISM OF HEARING AS REVEALED THROUGH 
EXPERIMENT ON THE MASKING EFFECT OF THERMAL 
NOISE 


By HARVEY FLETCHER 
BELL TELEPHONE LABORATORIES 


Read before the Academy April 26, 1938 


Considerable information concerning the dynamical properties of the 
hearing mechanism can be obtained from physical measurements on audi- 
tion. In fact, probably more precise information can be obtained by such 
measurements than from data obtained by animal experimentation, or 
from data on post mortem sections of human temporal bones. 

It is well known that when a person is immersed in a noise his ability 
to hear other sounds is decreased. The noise is said to produce a masking 
effect. Before describing the experiments on the masking effect of thermal 
noise I will first describe the method of obtaining thermal noise and how 
it is quantitatively defined. 

In an electrical conductor there is a statistical variation of the electrical 
potential difference between its two ends which is due to the thermal agita- 
tion of the atoms, including the electrons. This fluctuating voltage may 
be amplified by means of a vacuum tube amplifier. Such an electrical dis- 
turbance is called an electrical thermal noise. If this electrical thermal 
noise is sent through a telephone receiver, or a loud speaker, it is then con- 
verted into acoustical thermal noise. Such an acoustical thermal noise 
was used in the masking experiments discussed in this paper. 

Let I, be defined as the intensity per cycle of thermal noise using 10~'¢ 
watts* per square centimeter as a unit of intensity. If AJ is the intensity 
due to the thermal noise in the band of frequencies between f and f + Af, 
then 


AI 
L=—. 1 
(1) 
The spectrum level B in decibels is defined by 
B = 10 log J,. (2) 


In general J, and B will vary with the frequency. The curve showing the 
variation of B is called a spectrogram of the noise. 

Let us now examine more closely the relationship between the intensity 
of the noise and its masking effects. When a thermal noise having a given 
spectrogram is impressed upon the listener’s ears it causes the hearing 
mechanism to vibrate some parts more and other parts less. These vibra- 
tions stimulate the nerve endings. The ear has a selective action for 
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sounds of different frequencies so that the amount of agitation or stimulation 
given any set of nerve endings will depend upon their position. 

It is well known that the nerve endings responsible for auditory effects 
are located along the basilar membrane. According to the measurements 
of anatomists this membrane varies in width from 0.2 millimeter at the oval 
window end to 0.5 millimeter at the helicotrema end, and is about 30 milli- 
meters long. It will be seen that it is about 100 times as long as it is wide. 
Therefore, a single codrdinate x will define with sufficient accuracy the posi- 
tion of a patch of nerve endings. Then, if we take one per cent of the length 
it will contain a patch of nerve endings about as long as it is wide. 

To make this concept more concrete, a spiral line having approximately 
the same shape as the human cochlea along which the basilar membrane 
runs is shown in figure 1. The numbers along the spiral line give the ap- 


20 
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FIGURE 1 


23 24 


proximate distances in millimeters as you go from the helicotrema to the 
oval window. These dimensions vary in the ears of different persons. 
Also, the density of the nerve endings, that is, the number of endings per 
millimeter, is not accurately known. To avoid the necessity of using un- 
certain anatomical data, the position coérdinate x will not be taken in 
millimeters from the helicotrema but as the per cent of total nerve endings 
that is passed over in going from the position of maximum stimulation for 
the lowest audible frequencies (the helicotrema) to a position designated 
by the coérdinate x. The position coérdinate thus defined is equal to zero 
at the helicotrema end and equal to 100 at the oval window end. For ex- 
ample, as will be shown later, the position x of maximum stimulation for a 
pure tone of frequency 5000 cycles per second is equal approximately to 75. 
This means that 75 per cent of the nerve endings are in the direction toward 
the helicotrema and 25 per cent in the direction toward the oval window. 
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Let us designate the distance in millimeters from the helicotrema as 
x’. There is then a simple relation between the position coérdinate x and 
the distance in millimeters x’ if we know the density o, that is, the number 
of nerve endings per millimeter at each distance x’. This relation is given 
by 


dx" 


en 


MASKING LEVEL, Bm 


-10 10 20 30 40 50 
SPECTRUM LEVEL,B 
FIGURE 2 


In our discussions we will use the coérdinate x which can be determined 
directly from auditory experiments. 

When a continuous noise is impressed upon the ear it produces a certain 
stimulation at the position x. The stimulating agent will be taken as the 
vibrational power acting upon the nerve endings. Let the amount of this 
power between positions x and x + Ax be designated as AJ. Then the 
power J, for one per cent of the nerve endings is given by 
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J,= —- (4) 


The stimulation level S at any position x is defined by 
S, = 10 log Je (5) 
J, 

where J, is the threshold stimulation. Since here we are dealing with a 
ratio of vibrational powers, the uncertainty of determining whether to use 
amplitude velocity or acceleration or any function of them as the stimu- 
lating agent does not arise. The value S, can then be taken as the 
stimulation of the one per cent patch of nerve endings at the position x. 
This stimulating power can be spread over a small or a large patch of 
nerve endings and it will still produce just the threshold value. This is 
true because our experiments on loudness have shown that the nerve dis- 
charges near the threshold are proportional to the intensity of the sound. 

Let the following experiments be performed. A broad band of thermal 
noise whose spectrum level B is known at each frequency is impressed upon 
the ear of an observer. While this noise is present, a pure tone of known 
frequency f is raised to the intensity J,, so that the observers just perceive 
its presence. The intensity level of this tone which is masked will be desig- 
nated £,, and it is related to J,, by the equation 


= 10 log J,,. (6) 


For a fixed frequency a set of values of 8,, and B were obtained using differ- 
ent intensity levels of the noise. 

The experimental apparatus and methods have been described elsewhere. 
A typical set of results is shown in figure 2 for two frequencies 9500 cycles 
per second and 480 cycles per second. The spectrum level B of the noise 
at each of these frequencies is plotted as abscissae and the intensity level of 
the pure tone just masked by the noise is plotted as ordinates. It was found 
that this relation was independent of the frequency width of the band as 
long as this exceeded a critical value ranging from 100 c. p. s. in the low 
frequencies to 1000 c. p. s. in the high frequencies. Points above 6, = 
80 and below 6,, = 20 are usually above the straight line which goes 
through the points between these levels. There is a reason for both of these 
departures which cannot be dealt with here. It is the straight portions of 
the curves which give us the information we want to use in this discussion. 
Since they are straight lines, equations of the form 


6, = B+ K (7) 


will represent them, where K is a parameter varying with the frequency 
regions used. It is this variation that enables us to determine the relation 
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of the position of maximum stimulation x, to the frequency f of the im- 
pressed tone. Eleven frequency regions were tested and the results plotted 
on eleven straight lines like those of figure 2.+ The resulting values of K 
are shown in figure 2 for the various frequencies used. The relation be- 
tween the acoustic pattern and the corresponding stimulation pattern on 
the nerves is illustrated in figure 3. A 500-cycle band of noise between 500 
and 1000 is spread into a wide stimulation pattern of 16 per cent of the 
nerve endings, while this same band width between 7500 and 8000 of the 
same acoustic intensity acts upon only one per cent of the nerve endings at 
position 88 as indicated. Consequently, one must raise the level of a pure 
tone in this high frequency region to a higher value to be perceived than if 
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one were perceived in the low frequency region where the stimulation per 
nerve is not nearly so great. Conversely, if the masking level, compared to 
the spectrum level, is higher for the high frequency region than for the low 
frequency region, it shows that the frequencies are less spread out in the 
high and more spread out in the low frequency regions. For example, in 
the case illustrated the masking tone would be sixteen times greater for one 
than for the other. 

Let us consider this relation in a quantitative way. The intensity AJ 
in the thermal noise band between f and f + Af goes to stimulate mostly 
those nerve endings between x and x + Ax, where x corresponds to the 
position of maximum stimulation for a pure tone having a frequency f, 


Y 

| 
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and x + Ax to the position of maximum stimulation for a pure tone having 
a frequency of f+ Af. In other words, the band Af located at the frequency 
position f will be chosen so that the power in this particular noise band is 
spread over the particular patch of nerves Ax located at the position x. 
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As illustrated in figure 3 the patches adjacent to Ax will also be stimulated if 
a single band Af is impressed upon the ear. However, for any type of noise 
whose spectrum is not changing rapidly with frequency these adjacent 
patches will already be stimulated by the adjacent bands. For this type 
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of noise it is fairly accurate to assume that all the intensity in the band Af 
goes to stimulate the nerve patch Ax. The response of the nerves stimu- 
lated by small bands of thermal noise has been calculated. The re- 
sults of this calculation are shown in figure 4 for three different positions. 
It is seen that the ear is very selective with frequency and that the above 
assumption is a very good approximation. 

Since the thermal noise band and also the masked tone are in the same 
frequency region, they will have approximately the same frequency dis- 
crimination produced by the electro-acoustical system generating the sounds 
and also by the transmission system in theear. The acoustical intensity AJ 
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produces the stimulation AJ. The acoustical intensity J,, produces the 
stimulation J,, at the same nerve position x as for AJ. Therefore in the 
range where the system is linear 


(8) 


If we consider Af and Ax as differentials and integrate we have 


x= df. (9) 


It is now assumed that Jn is a constant C independent of position x and 
z 
consequently of the frequency f. 
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cC=—. (10) 


Also, equation (9) must fulfil the condition that x = 100 whenf = ~. If 
these conditions are applied and the values of B and £,, from equations (2), 
(6) and (7), then equation (9) reduces to 


—K/10 
fi df (11) 
0 
100 100 
12) 
Ly poms ( 


So we can substitute the values of K from figure 2 and integrate these equa- 


where 


tions graphically. 
In the region where the system is linear 
J, 
13 
(13) 


where J) is the threshold intensity of a band of thermal noise of small width 
located at the frequency position f corresponding to that of the tone which 
is masked, and Jp the stimulation necessary for threshold at the position x 
corresponding to f. If we multiply the numerator and denominator of the 
left-hand side of this equation by J, and then reduce the equation to deci- 
bels by taking 10 log of both sides, there results . 


Bm — Bo = S+ 10 log = C. (14) 
Zz 

The left-hand side of this equation is, by definition, the masking M. 
This shows that for a linear system the assumption that J,,/J, is a constant 
is equivalent to assuming that for a constant masking there is a constant 
stimulation. In other words, the amount of masking M is a measure of the 
stimulation S and the difference between these two quantities is simply 10 
log of the constant given by equation (10). If we express equation (7) in 
decibels and then rearrange the terms, there results 


By = B+ 100g C + 1010 2. (15) 
x 
From this equation it is seen that the intensity level 8,, of the masked tone 


can be calculated from the spectrum level B of the noise provided the rela- 
tion between x and f is known. 
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Using the values of K shown in figure 3, equations (11) and (12) were 
integrated graphically. The value for C was found to be 1.1, or within 
the experimental error equal to unity. This means that when the stimula- 
tion J,, due to the tone being masked is equal to the stimulation J, upon 
one per cent of the nerves then it will just be perceived by a typical observer. 
Also, this indicates that the masking M is directly equal to the stimulation 
S and may be considered a means of measuring it. The resulting values 
of f and x obtained from equations (11) and (12) are shown in the curve of 
figure 5. The relationship between f and x indicates that frequencies be- 
low about 150 cycles per second all have the position of maximum stimula- 
tion on the first one per cent patch of nerve endings. When the frequency 
of the stimulating tone reaches 1000 cycles per second the position of maxi- 


2000 


FIGURE 7 


mum stimulation has passed over approximately 30 per cent of the nerve 
endings. 

It will be remembered that the values of x refer to nerve patches in such 
a way that for equal intervals Ax there are an equal number of nerve end- 
ings. If the nerve endings are distributed uniformly along the basilar 
membrane then the values of x used in this discussion will correspond to 
one per cent of the total length of the basilar membrane. If, however, the 
nerve endings are not distributed uniformly, then the length dx must be 
stretched or contracted depending upon the nerve density being less or 
greater than the average. Let x’ be the distance from the helicotrema in 
millimeters. 

Guildt has made measurements of the number of ganglion cells along the 
basilar membrane. He divided the cochlea into four parts, namely (1) 
upper middle and apical, (2) lower middle, (3) upper basal, (4) lower basal. 
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He made a count of the total number of ganglion cells in each of these 
parts. They are given below in column 3 of table 1. Steinberg has esti- 


TABLE 1 

PosITIoNns x’ Count 
Upper Middle 

plus apical 0-7 4,200 0-14.5 
Lower Middle 7-12 5,600 14.5-33.8 
Upper Basal 12-19 8,500 33 .8-63.1 
Lower Basal 19-30 10,700 63. 1-100 

Total 0-30 29,000 0-100 


mated the distance in millimeters along the basilar membrane for each of 
these parts to be that given in the column under x’. The corresponding 
values of x are then easily calculated as the per cent of the total cells up to 
the position corresponding to x’. They are given in the last column. 

It is seen that these data on the number of ganglion cells allow us to calcu- 
late only three points besides the two end-points. These values of x and x’ 
are plotted on figure 6. A smooth curve passed through these points. 
Using this relation, then the values of x and f from figure 5 can be transferred 
to values of x’ and f. In figure 7 this relation is plotted on a spiral having 
the dimensions shown in table 1. The positions given in this spiral agree 
with previous determinations within the accuracy of those determinations. 


* All intensity values used in this paper are expressed by this same unit. 


+ I am indebted to Mr. W. A. Munson for doing this experimental work. 
t Guild, Stacy R., Acta Oto-Laryng., 17, 207-245 (1932). 


THE DIMENSIONLESS CONSTANTS OF PHYSICS 
By ARTHUR E. Haas 


UNIVERSITY OF NoTRE DAME 
Communicated June 8, 1938 


Physics appears to be characterized by several universal constants which 
are pure numbers. One is the ratio of the mass of the proton to the mass 
of the electron, i.e., m,/m = 1835, another the ratio e*/(h c), where e is the 
fundamental charge, h Planck’s constant and c the velocity of light. A 
third is the ratio e®/(G m, m), where G represents the gravitational con- 
stant. An important fourth number, under the assumption of a finite uni- 
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verse, is given by the ratio M/m,, where M denotes the total mass of the 
! universe. It appears that a fifth constant, the ratio of the total number of 
neutrons to the total number of protons, should be introduced, this ratio 
being referred either to a static universe, or under the assumption of an 
expanding universe, to its initial equilibrium state. 

It will be the task of this consideration to derive the result that the em- 
pirical constants of atomic physics yield for this ratio a value equal to unity 
to a high degree of approximation, provided that we can accept as correct 
the results of two earlier papers based on wave-mechanical deductions and 
published by Eddington and by Sitte and Glaser, respectively. 

Eddington,' in 1931, derived the simple relation 


R’= Na’, (1) 


where R denotes the equilibrium radius of the universe, N, the total num- 
ber of electrons and a (= e?/(mc?)) the classical radius of the electron. 

Sitte and Glaser? arrived at a formula which may be written in the follow- 
ing form 


h/(Mc) = (N, + N,)-, (2) 


where N, is the total number of protons; the existence of neutrons was not 
taken into account by Sitte and Glaser. Since, however, the only proper- 
ties of protons and electrons which play a réle in Sitte and Glaser’s consid- 
erations are that the mass of the electron is insignificant in comparison with 
the mass of the proton, and that the universe as a whole is electrically neu- 
tral, equation (2) may be generalized as follows: 


h/(Mc) = RN“ (3) 


where N is the total number of primordial particles of which N’ are heavy 
in the sense that the mass of the other (V — N’) particles may be disre- 
garded in comparison with their mass.* 

If we denote the number of neutrons by N, and remember that N, 
equals N,, we have 


N= N,+ 2N,, N’=N,+N, M = (N, + N,)m,. (4) 
If we put 
N, =x N, (5) 
we thus find 


h/(m, c) = (1 -+ 2)" (24 (6) 


= 
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By combining this relation with equation (1) and remembering the equal- 
ity of N, and N,, we obtain 


(2 + x) (1 +x)~'* = (m,/m) X e?/(he). (7) 


The empirical value of the right-hand side is 2.13 and the solution of equa- 
tion (7) is a value of x slightly above unity. It would just be equal to one 
if we could put the right-hand side of equation (7) equal to 2.12, thus di- 
minishing its empirical value by half a per cent. 

Thus, if the equations derived by Eddington and by Sitte and Glaser 
hold, we arrive at the conclusion that in the state of equilibrium the num- 
ber of neutrons nearly equals the number of protons. On the other hand, 
a relation, first pointed out by the author,‘ which connects purely atomic 
constants and which is fulfilled to within an error of about half a per cent, 
i.e., 


m p/m X e*/(hc) = 3//2 (8) 


may be derived from the formulae of Eddington and Sitte and Glaser by 
assuming the equality of the numbers of neutrons and protons in the 
equilibrium state. 

1A.§S. Eddington, Proc. Roy. Soc. London, A135, 605 (1931). 

2 K. Sitte and E. Glaser, Zeitschr. Physik, 88, 103 (1934). 


3 Concerning the details, cf. the paper by Sitte and Glaser. 
4A. E. Haas, Science, June (1938): The classical radius of the electron is related to 


the Compton wave-length of the proton as 3 to Vv 2. 


ON THE RELATIVE CONTRIBUTIONS OF NATURE AND 
NURTURE TO AVERAGE GROUP DIFFERENCES IN 
INTELLIGENCE 


By BARBARA S. BuRKS 
CARNEGIE INSTITUTION OF WASHINGTON, COLD SPRING HARBOR 


Communicated June 13, 1938 


It has become a more or less accepted commonplace that environmental 
differences, although having a demonstrable influence upon individual dif- 
ferences in IQ, are less potent than natural differences under the conditions 
ordinarily met in our urban culture. No wholly satisfactory quantitative 
estimate of the relative réles of nature and nurture has yet been possible, 
but various estimates agree in placing the contribution of nurture under 
50 per cent, and probably considerably under. (See references 1 to 5.)* 


i 
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Far less satisfactory is the present position of the problem of average 
socio-economic group differences in intelligence—in fact it is only in the last 
few years that it has been generally recognized and discussed as a separate 
problem requiring its own techniques for solution. There is no simple cor- 
respondence between the contributions of nature and nurture to group and 
individual differences, but the same types of data are crucial for both prob- 
lems. The writer has drawn upon two sources of data applicable to the 
problem of average group differences in IQ with respect to father’s occupa- 
tion—her previous study conducted at Stanford University of the IQ’s 
of foster children and ‘“‘own’’ children in relation to parental intelligence 
and home background, and the more recent study by Leahy at University 
of Minnesota dealing with the same type of material. 

The comparisons involved are unusually straightforward, for it can be 
shown that the difference in mean IQ of ‘‘own’”’ children grouped accord- 
ing to father’s occupation is composed of two additive heredity and environ- 
ment factors, the latter being given by the corresponding difference in 
group means of foster children. Wright has kindly furnished a proof, 
which depends upon the use of path coefficients, as follows: 


Child’s IQ can be represented as completely determined by the two 
factors—heredity (i.e., the child’s genetic constitution) and total environ- 
ment, which may be (and undoubtedly are) more or less correlated with 
each other. Occupational status of father is clearly correlated with the 
child’s total environment. The increased differences in the control data 
indicate that it is also correlated with the child’s heredity in the controls. 
The reasons (involving father’s intelligence) need not be represented. The 
following diagram represents these essential points: 


Occupation C 5 5 
Group Zotal Heredity 
Environment 


Teo = Re Veo Pew (control) 


tT 
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The observed differences in IQ in relation to given differences in occupa- 
tion group are of the nature of regression coefficients. 


beo = — = Pcutuo) — (control), 
= = (b' ce?’ zo) (foster). 
The partial regression coefficients measuring the concrete effects of H 
and E on C should be the same in foster and control data in spite of differ- 
ences in correlations, path coefficients or standard deviations. 


Per = P' ce Pou = 


If now we assume that the correlation between total environment and 
occupation group is the same in the two bodies of data but that there is no 
correlation between child’s heredity and occupation group in the foster 
data (no effective selection of children) : 


Tso = = 0, 
bco = ro— Pou’ (control), 
. 
where Pcetzo — = b'co given by the foster data. 


Thus the difference in the controls can properly be analyzed into two 
additive portions tracing through the correlation of occupation group with 
child’s environment, which is the same as the total difference in the foster 
data, and tracing through the correlation with child’s heredity. 


The published Minnesota study provides the needed tabulations for occu- 
pational group comparisons. It was necessary to go back to the original 
data of the Stanford study and tabulate them according to a similar 
scheme. 

The results are shown in tabular form. Table 1 contains the group 
means and dispersions of intelligence test scores; those of the parents as 
well as of the children are included by way of additional interest. Table 2 
contains derived estimates of the relative contributions of nature and 
nurture to occupational group differences. The “nurture” column of 
table 2 is simply the ratio (in per cent) of group difference in foster children 
to group difference in control ‘‘own’”’ children. The “‘nature’’ column is 100 
minus this ratio, or the ratio of increment of control difference to control 
difference. In table 3 are presented estimates of the relative contribu- 
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tions of nature and nurture based upon the total data of the Stanford and 
Minnesota studies, singly and in combination. Sampling errors of these 
estimates are also included. 

Two methods, giving virtually the same results, were used for obtaining 
combined estimates of the nature and nurture contributions (table 3). 

1. The occupational group differences shown in table 2 were summed 
for the control and foster groups separately, and the ratio of the foster 
sum to the control sum was computed. This procedure is sufficiently ac- 


TABLE 1 
MEANS AND DISPERSIONS OF INTELLIGENCE SCORES BY OCCUPATIONAL GROUP 


STANFORD STUDY MINNESOTA STUDY 


OCCUPATION OF FATHER OCCUPATION OF FATHER 
(OR FOSTER FATHER) (OR FOSTER FATHER) 
HIGHER SKILLED SLT. 
BUS., TRADES SKILLED 
PRO- SEMI- LOWER SKILLED BUS. & CLERI- SEMI- & DAY 
FESS. PROF. BUS. LABOR PROFESS. MGR. CAL SKILLED LABOR 
Il, Ill. 3¥. 1. II. Ill. Iv. v. 

Foster children 
Mean (IQ) 109.1 108.6 108.0 104.6 112.6 111.6 110.6 109.4 107.8 
S. D. (IQ) 17.2 14.8 16.7 
No. 32.0 47.0 41.0 48.0 43.0 38.0 44.0 45.0 24.0 
Control children 
Mean (IQ) 118.7 118.5 115.5 106.1 118.6 117.6 106.9 101.1 102.1 
S. D. (IQ) 15.4 12.2 18.6 12.4 ‘32:6 15:6. 14:3 21.0 
No. 18.0 338.0 27.0 18.0 40.0 423.0 48.0 46.0 23.0 
Foster parents* 
Mean 221.8 207.3 201.2 1841.7 58.6 59.6 49.6 39.7 38.4 
D 22.6 30.8 29.7 30.3 8.0 6.7 12.8 11.2 
No 24.0 40.0 34.0 34.0 
Control parents* 
Mean 221.6 221.8 192.0 176.2 64.6 657.1 51.8 44.0 38.3 
S. D. 24.4 30.4 33.4 31.6 §.4. 17:6 284 9.0 
No. 18.0. 32.0. 27:0 16.0 


* In the case of the Stanford study, data are for Stanford-Binet mental age in months 
of foster fathers and control fathers. In the case of the Minnesota study, data are for 
mid-foster parent and mid-parent point score on the Otis Test of Mental Ability. 


curate for an over-all appraisal when the population numbers and standard 
deviations of the sub-groups do not differ in an extreme manner. The 
sampling errors of the estimates thus obtained, however, are difficult to de- 
termine because of the lack of independence of the group differences en- 
tering the sums. 

2. Using the relation,{ 


Xp Ip Xo Zo 
oe 
e 
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and assuming a true r of unity (i.e., a constant proportional effect upon IQ 
of nurture and nature in the small range of IQ levels with which we are 
concerned), we may estimate the proportional contributions by the ratio** 


(2) 


The sampling variance of R is found by squaring and averaging logarith- 
mic differentials, viz., 


V,(R) = 


+ 2 


V (cp) (3) 


Op, Sc represent the dispersions of sets of scores in which each score is 


drawn from a unique universe of given central tendency and sampling 


TABLE 2 


ESTIMATES OF THE RELATIVE CONTRIBUTION OF NATURE AND NURTURE TO DIFFERENCES 
IN MEAN INTELLIGENCE SCORE OF CHILDREN GROUPED ACCORDING TO FATHER’S 


OccuPATION 
STANFORD STUDY MINNESOTA STUDY 
GROUPS CONTROL CONTRIB. CONTRIB. CONTROL CONTRIB. CONTRIB. 
COM- MINUS OF OF MINUS OF OF 


PARED CONTROL FOSTER FOSTER NURTURE NATURE CONTROL FOSTER FOSTER NURTURE NATURE 


I-III 3.2 | 2.1 34.4 65.6 11.7 2.0 9.7% 17.1. 82:9 
I-IV 12.6 4.5 8.1 35.8 64.2 17.5 32. 4:3 18:3 81.7 
16.5 4.8 11.7 29.1 70.9 
II-III 3.0 0.6 2.4 20.0 80.0 10.7 1.0 9.7 9.3 90.7 
II-IV 12.4 4.0 8.4 32.3 67.7 16.5 2.2 14.3 13.3 86.7 
II-V 15.5 3.8 11.7 24.5 75.5 
III-IV 9.4 3.4 6.0 36.2 63.8 5.8 1.2 20 204-80 
4.8 2.8 2.0 58.4 41.6 


fluctuation. Hence their sampling variances (to be used in formula 3) 
cannot be calculated by the usual formula. A formula for V,(c,), V,(o¢) 
may be derived as follows: 


(x, — 2)? 

(4) 
(x, — #)? 
n—1 


Taking differentials, squaring and averaging, 


Vo) = 5 Vole) + 


— 1)? 


[V.(x,) = o*x,/n; obtainable from table 1]. 
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The difference in estimated contributions of nature and nurture yielded 
by the Stanford study (<2/3,>1/3) and the Minnesota study (>3/4, 
<1/4) may be due to sampling error alone, although examination of the 
slope of the regression of foster means upon control means in the two studies 
suggests that the estimate of the potency of nurture is rather consistently 
lower in the Minnesota study than in the Stanford study. This may per- 
haps be accounted for by the fact that there was less opportunity for ‘‘se- 
lective placement” in the Minnesota foster children, who were limited to 
those placed in adoptive homes under 6 months of age, than in the Cali- 
fornia foster children who were placed up to 12 months of age. 


TABLE 3 


ESTIMATES OF THE RELATIVE CONTRIBUTION OF NATURE AND NURTURE BASED UPON 
Tora DATA, WITH TESTS OF SIGNIFICANCE 


STAN- MINNE- 
FORD STUDY SOTA STUDY 
Unweighted sum of control group differences 40.8 99.0 
Unweighted sum of foster group differences 14.1 23.6 
. diff. 
Estimated contribution of nurture 0.345 0.238 
control grp. diff. 
Estimated contribution of nature (1-nurture contrib.) 0.655 0.762 
Dispersion of control group means (¢¢) 5.92 8.37 
Dispersion of foster group means (cp) 2.03 1.87 
Estimated contribution of nurture (R = of/oc¢) 0.344 0.224 
Estimated contribution of nature (1 — R) 0.656 0.776 
Sampling variance of R (Vo(R)) 0.054 0.017 
Sampling error of R (V V.(R)) 0.232 0.131 
Information (1/Vo(R)) 18.59 58.82 
Combined R (weighted according to information) 0.253 
1-R 0.747 
V.(R) 0.013 
Sampling error (VV Vo(R) 0.114 


Combining the estimates in the two studies according to the amount of 
information yielded by each (inversely as their sampling variances), we 
arrive at an estimated 3/4, 1/4 as the relative contributions of nature and 
nurture, with approximately even chances that the contribution of nurture 
is in truth between 18 and 33 per cent. 

Finally, it may be pointed out that differences in mean intelligence scores 
of the parents according to occupational group might be expected to derive 
from natural differences in a higher proportion than do those of their off- 
spring, since intelligence is one of the components that enters directly into 
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the selection of an occupation, or the selection of a person by an occupation. 
This truism is borne out by the higher regression of intelligence on occupa- 
tion in the parent groups than in the offspring groups. 

* Reference 3 contains quantitative estimates of the contributions of nature and nur- 
ture which have been strongly criticized. However, the simple correlation, 0.67, be- 
tween the intelligence scores of identical twins reared apart—which would rise even 
higher if corrected for attenuation, and for the restricted range of talent—is in itself 
weighty evidence for the predominant influence of natural inheritance. 

t To avoid the slightest possibility of bias, the classification of fathers’ occupations 
from the Stanford data was made “‘blind,” i.e., with no knowledge as to the intelligence 
of the fathers, mothers, their children or foster children. 

} In the formulae which follow, x, x and o refer to the occupational group means, the 
mean of the means and the standard deviation of the means, respectively. 


** This formula neglects a small error factor involving the magnitude of the obtained 
o’s in relation to the true o’s [i.e., ¢,,2 = o°(r), where 7 is a reliability coefficient ]. 

1 Burks, B. S., 27th Yearbook Nat. Soc. Stud. Educ., Part I, 219-316 (1928). 

2 Leahy, A. M., Genet. Psychol. Monog., 17, 235-308 (19385). 

3 Newman, H. H., Freeman, F. N., and Holzinger, K. J., Twins, Univ. Chicago Press, 
pp. xvi + 369 (1937). 

4 Willoughby, R. R., 27th Yearbook Nat. Soc. Stud. Educ., Part I, 55-59 (1928). 

5 Wright, S., Jour. Amer. Stat. Assoc. Suppl., 26, 155-163 (1931). 


A METAGALACTIC DENSITY GRADIENT 
By HARLOW SHAPLEY 
HARVARD COLLEGE OBSERVATORY 


Communicated June 14, 1938 


1. Ina discussion of the evidence bearing on the isotropic distribution 
of galaxies it was shown four years ago that at a distance of about 108 
light-years to the north of the galactic plane the density is thirty or more 
per cent greater than at the same distance to the south.’ Several different 
ways of discussing the extensive observational material yielded the higher 
density in the north. Similar inequalities were found to be present also at 
shorter distances, but beyond 10* light-years both the Harvard and the 
Mount Wilson data suggest tentatively the disappearance of the north- 
south inequality. 

Quantitatively the north-south difference is not closely determined by 
the early work because the nebular counts were made in relatively small 
sample areas; but the existence of the difference is unquestionable, as is 
also the occurrence of non-uniformity in nebular distribution arising from 
the metagalactic clouds.? These well-established large-scale inequalities 
have been pretty generally ignored in discussions of the expanding uni- 
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verse; the assumptions of isotropy and large-scale uniformity have per- 
sisted, leading of late to somewhat remarkable deductions. 

From surveys with reflecting telescopes Hubble has found a radial density 
gradient which, over an interval of 2.5 X 108 light-years, amounts to a den- 
sity change of something less than twenty per cent.* To eradicate this ap- 
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FIGURE 1 
The density distribution in a thirty-degree zone across the south galactic pole. 


parent increase in density with distance and restore nebular distribution 
to essential uniformity, he suggests the abandonment of the interpreta- 
tion of red-shift as velocity of recession and the abandonment, therefore, 
of the hypothesis of the expanding universe, introducing in its place some 
new principle to account for the observed red-shift. But the evidence for 
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a radial gradient is not very strong. It can be readily shown that if the 
Mount Wilson surveys are analyzed separately for the two galactic hemi- 
spheres the gradient to the north is inappreciable. 

In a recent discussion of the distribution of over seventy-five thousand 
galaxies in the south galactic polar cap‘ we have shown, in another way, 
that the accurate determination of a radial gradient is not yet possible be- 
cause of the large density irregularities that exist up to distances of 108 light- 
years, if not throughout all attainable space. If, combining the Mount 
Wilson and Harvard counts of galaxies, we use as a standard for our density 
calculation the mean value obtained from the census in the first (northwest) 
longitude quadrant of the polar cap, we find that the density increases 


2.00 


1.80 


1.00 


B—30° -50° —70° Pole -70° -—50° -—30° 
FIGURE 2 
The metagalactic gradient as shown by a survey using the 
nine central square degrees of one hundred plates. Numbers 
of galaxies have been reduced to a common magnitude limit. 


conspicuously with distance; if we use as standard the comparable material 
from the third quadrant, the density decreases slightly. These quadrants 
contain over 39,000 observed galaxies, and the determination of the differ- 
ence in the density therefore has high weight. 

2. The survey of the south galactic cap provides evidence of a long-dis- 
tance transverse gradient or inequality in metagalactic density* that is as 
pronounced as the north-south inequality referred to above. It is more 
securely established by the observations because of the greater amount 
of material involved. Since this progressive inequality observed in the 
south galactic cap is across the line of sight, its measurement is free of first 
order difficulties with magnitude scales. The radial gradient suggested by 
the Mount Wilson surveys directly assumes a correct magnitude scale, and 
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disappears, according to Hubble,’ if the scale is too restricted by 0.12 be- 
tween magnitudes 18.5 and 21.0. The north-south inequalities found from 
the Harvard surveys four years ago involve the inter-comparison of fre- 
quencies in opposite regions, and such nebular counts are also moderately 
free of the effects of scale errors in the magnitude system. 

3. The work on the south galactic cap (galactic latitudes —55° to —90°) 
has now been extended in order to explore further the progressive though 
somewhat irregular increase of nebular frequency as the survey crosses the 
south galactic pole from the northwest to the southeast. The increase was 
grossly shown when we found in the four quadrants the following numbers 
of galaxies,* to the average magnitude 18.2 (number of plates in paren- 
theses) : 


Northwest 10,800 (18) Southwest 22,800 (21) 
Northeast 19,000 (18) Southeast 23,000 (23) 


Reducing the counts to a common magnitude and using only results from 
the inner nine square degrees on each plate, we obtain the following values 
for the number per square degree: 


Northwest 31.87 Southwest 51.10 
Northeast 54.17 Southeast 66.54 


The richest quadrant, on the average, is twice as rich as the poorest; the 
transverse distance involved is of the order of 108 light years. 

4. The manner in which the survey has been extended is most clearly 
indicated by the polar diagram in figure 1. The south galactic pole is at 
the center; galactic latitudes and longitudes form the codrdinate system. 
In a zone thirty degrees wide across the south galactic polar cap, centrally 
aligned on the longitude circle 45°, 225°, the positions are plotted of one 
hundred survey plates (each of three hours’ exposure with the Bruce tele- 
scope). The frequency of galaxies is roughly indicated by the shading. 
The area of each circle is approximately twenty-five square degrees— 
the same as that covered effectively on the sky by a Bruce plate. 

The diagram shows that the density continues to increase to the south- 
east, beyond the limits of the earlier survey. Whether the maximum has 
been reached at latitude —35° (right ascension 5", declination —60°), be- 
fore the absorption of the Milky Way system makes itself appreciably felt, 
cannot yet be determined. The observed falling-off in density beyond the 
maximum is probably a feature of metagalactic structure in this region, but 
may be also in part caused by outlying Milky Way absorption, which has 

* The objects lying outside the twenty-five square degrees on each of the eighty plates 
are not included. They approximately balance in number the objects appearing twice 


in overlapping of plates for which correction has not been made in the summation 
above. 


286 ASTRONOMY: H. SHAPLEY Proc. N. A. S. 


been shown in a study of the south equatorial polar cap to be widespread in 
this part of the sky.” 

In table 1 mean values, reduced to the convenient magnitude limit 17.9, 
are given for the results from the central nine square degrees of the plates, 
and also for the surrounding sixteen square degrees. Each mean includes 
the results of five plates, and therefore forty-five square degrees and eighty 
square degrees, respectively, are represented by each entry and by the 
plotted points of figures 2 and 3. The last column of the table gives for 
each group of five plates the total number of galaxies actually observed 
within the twenty-five square degrees. 


TABLE 1 


PosITIONS AND AVERAGE NUMBERS OF GALAXIES PER SQUARE DEGREE 


GALACTIC _ LOGARITHMS _ TOTAL NUMBER 
LONG. LAT. Nas Nso OF GALAXIES 
45° —29° 1.14 0.88 3,630 

—36 1.10 0.85 4,353 
—47 1.18 1.00 3,788 
—61 0.99 0.87 2,773 
—68.5 1.07 0.92 3,475 
—76.5 1.12 0.90 3,440 
—81l 1.13 0.91 5,992 
—87 1.21 1.08 4,789 
225 —88 1.08 0.88 4,570 
—81 1.20 0.95 4,648 
—75.5 1.30 1.10 5,650 
—71 1.30 1.08 5,846 
—68.5 1.50 1.29 7,582 
—58 1.45 1.35 4,363 
—54 1.55 1.38 4,820 
—50.5 1.60 1.47 4,918 
—45.5 1.76 1.58 6,669 
—41.5 1.90 Beg | 6,557 
—35 1.91 1.70 5,706 
—28 1.66 1.48 3,978 
Total 97,547 


So much area and so many objects are included in each mean that they 
no doubt smooth out large local irregularities. But the plots are given here 
chiefly to show the large-scale non-uniformity that stretches over 125 de- 
grees of the southern sky. The completion of the eighteenth magnitude 
survey for the southern hemisphere will permit more detailed description 
of the nebular distribution; but as it stands, we see a conspicuous and 
steady rise in the frequency from Ns = 10, at about latitude —60° in the 
northwest quadrant, to Nas = 80, at about latitude —35° in the southeast 
quadrant. 
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The distance of the objects that chiefly establish the mean values is ap- 
proximately 1.3 X 10° light-years. From the region of minimum fre- 
quency to that of maximum frequency is therefore about 2 X 108 light-years 
—a transverse distance comparable to that over which Hubble found an 
apparent radial gradient. Since we find here the frequency increasing by 
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0.80) 
B—30° -—50° —70° Pole —70° -—50° —30° 
FIGURE 3 
The metagalactic gradient as shown by a survey using the 
sixteen square degrees surrounding the central 9-square on a 


m 
hundred plates. Galaxies only to a limit 0.3 brighter than 
that for figure 2 are represented by this density curve. 


several hundred per cent, it would appear reasonable to accept a twenty 
per cent radial density change as a minor detail of metagalactic structure, 
probably not associated significantly with the interpretation of the red- 
shift or with space curvature. 


1 Harv. Bull., 894 (1934). 

2 Harv. Ann., 105, Tercentenary Paper No. 8 (1937); see also Harv. Ann., 88, No. 5 
(1935). 

3 Mon. Not. R. A. S., 97, 508 (1937). 

4 These PROCEEDINGS, 24, 148 (1938); Harv. Circ., 423 (1937). 


oo 


288 PALEONTOLOGY: C. STOCK Proc. N. A. S. 


A TARSIID PRIMATE AND A MIXODECTID FROM THE POWAY 
EOCENE, CALIFORNIA 


By CHESTER STOCK 


BaLcH GRADUATE SCHOOL OF THE GEOLOGICAL SCIENCES, CALIFORNIA INSTITUTE OF 
TECHNOLOGY 


Communicated June 7, 1938 


Introduction.—Following the discovery and description of a titanothere! 
in the Poway sands of San Diego County, continued excavation in this im- 
portant Eocene horizon has produced small lots of mammalian remains. 
One of these included specimens of a tarsiid primate and a mixodectid. 
Worthy of mention is the fact that the former represents a third member of 
the Tarsiidae to be discovered in the early Tertiary of the Pacific Coast 
region. Since the Poway fauna possesses considerable significance in any 
attempt to correlate the later Eocene horizons of the Pacific Coast with 
those of the Cordilleran Province, it seems desirable to place on record 
these new mammals from the Poway. 


Yumanius woodringi, n. gen. and n. sp. 


Type Specimen.—Right maxillary fragment with M1 and M2, No. 2233, 
Calif. Inst. Tech. Vert. Pale. Coll., Plate 1, figure 1. 

Paratypes.—Fragment of right ramus with M1-—M3, No. 2234, Plate 1, 
figures 2, 2a and ramal fragment with M2, M3, No. 2235, Plate 1, figures 3, 
3a. 

Locality —White sandstones associated with the Poway conglomerate 
and exposed on west bank of San Diego River, approximately one-quarter 
mile north and east of San Diego Mission; C. I. T. Vert. Pale., Loc. 249. 

Generic and Specific Characters—Upper molars differ from those of 
Euryacodon lepidus in possessing well developed and clearly defined hypo- 
cone or postero-internal cusp. Intermediate cusps also well developed. 
A second small cusp, comparable in size to that of the metaconule is situ- 
ated between the latter and the base of the protocone. M2 smaller than 
that in Dyseolemur sylvestris and hypocone, in contrast to that in the latter, 
more like a cingular cusp. 

Lower molars 2 and 3 with trigonid portion of crown compressed antero- 
posteriorly and paraconid considerably reduced in size, more so than in 
Euryacodon. M3 with broader posterior rim to heel. Distinctly different 
from Dyseolemur in reduction in size of trigonids in posterior molars and in 
absence of metastylid. I take pledsure in naming the species for my 
former colleague, Dr. Wendell P. Woodring of the U. S. Geological Survey. 

Comparisons.—The upper molar teeth of Yumanius woodringi exhibit a 
cingulum on the outer, anterior and posterior sides of the crowns but the 
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basal ledge is lacking on the inner side. The small cuspule shown by 
Wortman? as situated at the inner base of the protocone in M2 of Euryaco- 
don lepidus is absent in the comparable tooth of Y. woodringi, but a slight 
development of the cingulum does occur in the second molar at the antero- 
internal base of the cusp. In the type specimen of E. lepidus, No. 11813 
Y. P. M., the enamel of the crowns of the upper molars is smooth, the 
postero-internal cusp is not clearly indicated and the intermediate cus- 
pules are small. In the molars of Yumanius, the postero-internal cusp 
shows pronounced development and the structural features of the crown 
become more complicated by the addition of a small cusp between meta- 
conule and protocone. In addition the anterior ridge of the protocone, di- 
rected toward the protoconule, is more strongly formed in Yumanius than 
in Euryacodon. 

The molars are unfortunately the only teeth preserved of the lower denti- 
tion. These resemble most closely the comparable teeth in Euryacodon. 
The enamel of the crowns is wrinkled, not smooth as in Omomys; the 
teeth are small, the paraconid while of normal size in M1 becomes reduced 
to vestigial proportions or to the point of disappearance and the trigonid 
parts of the teeth are compressed anteroposteriorly. This compression of 
the trigonid portion of the crown is greater in the Poway genus than in 
Euryacodon or Anaptomorphus. In M3 the posterior rim of the heel is 
broad and the end of the crown is not pointed as in Euryacodon. While 
no distinct metastylid is present as in Dyseolemur and Washakius, the pos- 
teriorly directed wing or ridge of the metaconid is extended in fore and aft 
line and is compressed transversely in Yumanius. 

In the type of Anaptomorphus aemulus Cope, No. 5010 A. M. N. H., the 
teeth are smaller than in Yumanius. In M1 of both forms the paraconid is 
well developed. This cusp becomes vestigial in M2, but its reduction and 
the anteroposterior compression of the trigonid is greater in Yumanius. 
In the lower molars of the Bridger genus the external cingulum is not so 
well developed as in the San Diego specimens. 

Discussion.—The upper and lower jaw specimens with teeth herein re- 
ferred to the new genus Yumanius were found at one locality in the Poway 
but were not directly associated. The view that the type and paratypes 
belong to two distinct tarsiid genera cannot be ignored, although apparent 
absence of evidence other than that of fortuitous occurrence makes this 
possibility a seemingly remote one. Presence of this material in the Po- 
way gives strength to the belief that Wortman’s reference of the three 
lower jaw fragments of a tarsiid from the Bridger to Euryacodon lepidus is 
correct. It follows, as Matthew* pointed out, that Anaptomorphus and 
Euryacodon are probably identical. Among the several tarsiid genera re- 
corded from the Bridger middle Eocene the Poway type is most closely 
related to the Euryacodon-Anaptomorphus stock. Yumanius appears to 
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have carried the line a stage farther and its presence at the San Diego 
locality suggests at least a post-Bridger age for the Poway. Curiously, Y. 
woodringi is more closely related to the Bridger species E. lepidus than to 
the Sespe type, Dyseolemur pacificus. 


Microsyops kratos, n. sp. 


Type Specimen.—A left ramus with Dp3, P4—M3 inclusive and root 
fragments of the anterior premolars and canine, No. 2232, C. I. T. Vert. 
Pale. Coll., Plate 1, figures 4, 4a. 

Locality—Poway sandstones and conglomerates; C. I. T. Vert. Pale., 
Locality 249. 

Specific Characters.—Differing from Microsyops elegans and M. annectens 
most noticeably in size. Smaller than Craseops sylvestris, but with lower 
molars, relative to size of jaw, larger than in latter species. 

Description and Comparison.—This is the largest species of Microsyops 
so far recorded, being distinctly larger than M. elegans and M. annectens 
and resembling in size the species of Craseops recorded from the upper 
Eocene stage of the Sespe of Southern California. Dp3 is a two-rooted tooth 
with simple crown. The latter is triangular in cross-section with apex 
placed forward. A downwardly directed ridge on the posterior surface of 
the principal cusp divides this surface into two parts. When the specimen 
was collected this tooth was loosely attached and on removal disclosed a 
small bit of enamel of a permanent tooth beneath. In front of Dp3 is the 
exposed cross-section of the anterior premolar and presumably the canine. 
P4 is the largest tooth of the series, P4-M3 inclusive. In this tooth only 
a vestige of the paraconid ridge remains and the metaconid does not reach 
the height of the protoconid. The basin of the talonid is deeply exca- 
vated. This tooth and the molars which follow all show a well defined 
cingulum at the base of the protoconid and below the notch between this 
cusp and the hypoconid. In M3 the posterior wall of the hypoconulid is 
broken away. The ramus is sturdy and of approximately same depth 


PLATE 1 


Yumanius woodringi, n. gen. and n. sp. 

Figure 1, type specimen, No. 2233, skull fragment with M2 and M3, occlusal view. 
Figures 2, 2a, No. 2234, jaw fragment with M1-M3, view of inner side and occlusal 
view. Note: In figure 2a, the paraconid of M1 has been restored. Figure 3, 3a, No. 
2235, outer and occlusal views. All figures approximately x6. 


Microsyops kratos, n. sp. 


Figures 4, 4a, type specimen, No. 2232, ramus with cheek teeth, lateral and occlusal 
views, X2. 


Calif. Inst. Tech. Vert. Pale. Coll. Poway Eocene, San Diego Co., Calif, 
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throughout its length. Two mental foramina are situated below the 
anterior premolars. 

Microsyops elegans from the Bridger, as represented by No. 12590, 
A. M. N. H., is a distinctly smaller, more slender form than the type from 
the Poway. In the former specimen the lower end of the symphyseal 
contact surface extends farther back than in No. 2232 from the Poway 
Eocene, and the second premolar is not crowded between adjacent teeth 
as in the latter. Moreover, there appears to be less discrepancy in length 


COMPARATIVE MEASUREMENTS (IN MILLIMETERS) 


Microsyops Microsyops Microsyops Craseops 
kratos elegans annectens sylvestris 
TYPE SPECIMEN TYPE SPECIMEN 


No. 2232 No. 12590 No. 11791 No. 1399 
A A. M.N. H. Y. P. M. 
Length from anterior end D#3 to pos- 


terior end of M3................ 928.7 
Length from anterior end P2 to pos- 

terior end of M3................ ie 20.4 
Length from anterior end P3 to pos- 

terior end of M3................ eee 18 


Length from anterior end P4 to pos- 
terior end of M3................ 


— 

i=] 

PASH WS aN 


Depth of ramusat posterior end of M3 

Depth of ramus at anterior end of M2 

Thickness of ramus below M2....... 

Thickness of ramus below M3....... 
° Approximate 
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4.5 


between P4 and M1 in the Bridger species than in that from the San Diego 
locality. 

No. 12050 A. M. N. H., also from the Bridger Eocene and determined as 
Microsyops annectens, approaches the Poway specimen in size. While the 
ramus exhibits a depth like that in the Californian type, its thickness is dis- 
tinctly less. 

No. 1399 C. I. T., representing Craseops sylvestris from the Sespe pos- 
sesses a more massive ramus with the upper anterior border and the anterior 
end of the masseteric area decidedly more pronounced than in Microsyops 
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kratos. The two molar teeth preserved in No. 1399, M2 and M3, are 
slightly broader in comparison to their length in Craseops than in Micro- 
syops from the Poway. . 


1 Stock, C., Proc. Nat. Acad. Sci., 23, 48-53 (1937). 
2 Wortman, J. L., Amer. Jour. Sct., Ser. 4, 17, 139-140, fig. 133 (1904). 
3 Matthew, W. D., and Granger, W., Bull. Amer. Mus. Nat. Hist., 34, 457 (1915). 


LARGEST DEGREE OF A SUBSTITUTION IN THE GROUPS OF A 
GIVEN DEGREE 


By G. A. MILLER 


DEPARTMENT OF MATHEMATICS, UNIVERSITY OF ILLINOIS 
Communicated May 28, 1938 


Let G represent a substitution group of degree and let s represent a 
substitution of G which has the property that none of the substitutions of G 
has a degree which exceeds that of s. It follows directly from the average 
number of letters in the substitutions of a given group that the degree of s 
cannot be less than n/2 + 1. It has recently been proved! that when the 
degree of s has this smallest possible value, then all the substitutions of G 
besides the identity are of the same degree and G is an abelian group whose 
order is a power of 2. Moreover, m must then be of the form 2” + 2” ~ }+ 
...+ 2 and there is one and only one such group for an arbitrary positive 
integral value of m. The order of this group is 2” and it is conformal with 
the abelian group of type 1”. Every substitution besides the identity of an 
arbitrary one of this infinite system of groups has exactly half of its letters 
in common with every other substitution besides the identity of this group 
whenever m > 2. 

When the degree of s exceeds n/2 + 1 it is at least as large as n/2 + 
3/2, and when it has this value there are only three possible groups, viz., 
the two groups of degree 3 and the nonabelian group of degree 5 and of 
order 6 whose transitive constituents are of degrees 2 and 3, respectively. 
As this fact was also established in the article to which we referred above, 
the smallest possible relative degree of s which has not been considered is 
n/2+ 2. When the degree of s is n/2 + 1 the group of degree 2 is the only 
possible transitive group, and when it is n/2 + 3/2 the two groups of 
degree 3 are the only possible transitive groups. It is obvious that when 
the degree of s is n/2 + 2 then the five transitive groups of degree 4 are the 
only possible transitive groups but there are then two intransitive groups 
in addition to three infinite systems of such groups whose degrees are of the 
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form 2" + 2”— 14 ...+.4,m > 2, in two cases, and of the form 2” + !in 
the third case. In one of the first two cases the order of G is 2” ~ ' while in 
the other itis 2”. In the third case the order of G is also 2”. 

Each of the possible groups of degree 4 has the property that its sub- 
stitutions of highest degrees are of degree /2 + 2. Hence it may be 
assumed in what follows that ” > 4 unless the contrary is stated. If at 
least one of the transitive constituents of G is of odd degree and the sub- 
stitution of highest degree in G is of degree n/2 + 2, then two of the transi- 
tive constituents of G are of degree 3. Each of the remaining transitive 
consituents of G is then of degree 2 and hence the average number of letters 
in the substitutions of G is n/2 + 1. The value of m cannot exceed 8 be- 
cause not more than four of the letters of the substitution of largest 
degree could appear in the group whose transitive constituents are 3, 3, 
and when » = & there are obviously two such groups of orders 6 and 18, 
respectively. Hence there are two and only two groups of degree m in 
which a substitution of highest degree is of degree n/2 + 2 and which in- 
volve separately at least one transitive constituent of odd degree. The 
transitive constituents of these groups are of degrees 2, 3, 3, respectively. 

It remains to determine the groups of degree ” whose substitutions of 
largest degrees involve n/2 + 2 letters and all of whose transitive constitu- 
ents are of even degrees. The largest value of such a degree is 4 and not 
more than one of the transitive constituents of G can be of thisdegree. If G 
contains one such constituent and » > 4 its remaining constituents are of 
degree 2 and hence s involves at least (n — 4)/2 + 1 letters from this con- 
stituent. This is impossible since it involves 4 letters from the transitive 
constituent of degree 4. It results therefore that when the degree of s is 
n/2 + 2 and each of the transitive constituents of G is of even degree then 
all of these degrees are equal to 2 whenever » > 4. It remains only to 
consider the cases when each of the transitive constituents of G is of degree 
2 and when G is therefore an abelian group. 

The fundamental infinite system of such groups can then be constructed 
by choosing for s a substitution of degree 2”, m > 2, and extending it by a 
substitution of the same degree which has half its letters in common with s 
just as in the case when the substitution of largest degree in G is of degree 
n/2+ 1. When m > 2 we extend this group of order 4 by a substitution 
which has half of its letters in common with each of the three substitutions 
of order 2 in this group of order 4. This process is repeated until each of the 
substitutions has four and only four letters in common with every other 
one while the separate pairs have half of their letters in common. The 
resulting group is of degree 2” + 2”~ ...+ 4 and of order 2” ~ 
All of its substitutions, besides the identity, are of degree 2”. By ex- 
tending this group by a substitution of degree 2” — 2 which has two letters 
in common with each of its substitutions but involves no additional letter 
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there results a group of order 2” and of degree 2” + 2"~ 14 ...4 4 
whose substitution of largest degree involves n/2 + 2 letters. 

Half of the substitutions of the latter group are negative and are of the 
same degree. The smallest two groups which belong to this system are the 
intransitive group of degree 4 and of order 4, and the group of degree 12 
and of order 8 which involves four substitutions of degree 6 but no sub- 
stitution of a larger degree than 8. To obtain the third infinite system of 
groups of degree n whose substitution of largest degree is of degree n/2 + 
2 we adjoin to an arbitrary group of the former of the two infinite sys- 
tems noted in the preceding paragraph a substitution of degree 2” + 4 
which has two letters in common with each of its substitutions but involves 
four additional letters. The smallest group which belongs to this system is 
of order 4 and of degree 8. There are therefore two infinite systems of groups 
of degree n in which the substitution of largest degree is of degree n/2 + 2 
whenever n is of the form 2™ + 2™~1+4 ...44. When n is of the form 
2™ +! there is an additional such system of groups. 

To prove that these three infinite systems of groups include all of the 
groups of degree in which the substitutions of largest degree are of degree 
n/2 + 2 and each of the systems of intransitivity is of degree 2, it is de- 
sirable to consider separately the three cases in which G is composed of the 
substitutions of highest degree together with the identity, in which the 
substitutions of highest degree generate a proper subgroup of G, and in 
which these substitutions generate G but G involves also substitutions which 
are not of highest degree. In each of these cases G involves at least n/4 
substitutions of highest degree since the average number of the letters in its 
substitutions is »/2 and hence it must involve 1/4 substitutions of highest 
degree to account for the fact that the identity is of zero degree, since each 
substitution of highest degree involves two more letters than the average. 

When all the substitutions of G besides the identity are of degree n/2 + 
2 then this degree is 2” since every two substitutions of order 2 contained 
in G have half of their letters in common. It therefore results that n = 
2” 4+2"-14 ...+4+ 4 and the order of G is 2”~*. Hence G is one of 
the groups of the first of the three infinite systems noted above whenever 
it is composed of its substitutions of highest degree together with the iden- 
tity. It is obvious that a substitution s which has two letters in common 
with each of the substitutions of order 2 in such a G is negative and of 
degree n/2. The products of s into all the substitutions of G give all of 
the substitutions which have this property with respect to a given G and 
hence there results a group which is completely determined by the given 
G and of twice its order. This group belongs to the second of the infinite 
system noted above and all of the groups which can be formed in this way 
constitute this system. 

It remains to consider the case when the substitutions of highest degree 
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in G generate G but generate also substitutions of lower degree. In this 
case the substitutions of highest degree are negative since this degree is a 
power of 2 increased by 2. In fact, if the degree of these substitutions 
is a power of 2 increased by 2k and we construct G as noted above the substi- 
tution of highest degree in the resulting group is of the form n/2 + k + 1. 
In this construction we extend successively the group formed by the sub- 
stitution built on the part of a substitution of largest degree whose degree 
is a power of 2. It therefore results that the three infinite systems of 
groups noted above include all the groups of degree » which involve only 
systems of intransitivity of degree 2 and whose substitutions of largest 
degree involve exactly m/2 + 2 letters. The orders of the groups of the 
last system are one-half of their degrees. 

When the largest degree of a substitution in a group of degree n exceeds 
n/2 + 2 it is at least as large as n/2 + 5/2 and when it has this value m is 
odd and the number of possible groups is infinite since if there is a transitive 
constituent of degree 3 there is an infinite number of groups on the remain- 
ing n — 3 letters whose substitutions of largest degree are (n — 3)/2 + 1. 
Hence there is an infinite number of groups on letters whose substitutions 
of largest degree are of degree (n — 3)/2 + 4 or nm/2 + 5/2. Such 
groups can be constructed by forming the direct products or by establishing 
an isomorphism between the symmetric group of degree 3 and a group 
whose order is a power of 2 and whose substitutions of largest degree are 
of degree (n — 3)/2 + 1. They can also be constructed by forming the 
direct products of the alternating group of degree 3 and the latter groups 
as well as by the direct products of the three groups of degree whose sub. 
stitution of largest degree involves n/2 + 3/2 letters. 

It has been noted that the smallest relative number of letters in a group of 
degree n appears only when 7 is of the special form 2” + 2”~ '+... +2, 
and that when 1 is of the form 2” + ' then there are groups whose substi- 
tutions of largest degree do not involve more than n/2 + 2 letters. From 
the latter of these theorems it results directly that for an arbitrary value 
of there is a substitution group of degree m which does not involve more 
than 3n/4 + 2 letters since there is a number which is a power of 2 and ex- 
ceeds n/2. This degree is usually much smaller as follows from repeated 
applications of the given results. 


1 Proc. Nat. Acad. Sci., 24, 202-204 (1938). 
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REMARKS ON RIEMANN’S DOCTORAL DISSERTATION 
By JEssE DouGLas 
_ BROOKLYN, N. Y. 
Communicated June 10, 1938 
1. The author has based the solution of Plateau’s problem on the mini- 


mum principle 


D(H) = (E + G)du dv 


among all harmonic vectors H(u, v) which represent a surface bounded by 
a given contour I in m-dimensional euclidean space.! The (u, v) domain 
may be supposed to be the unit circle C. 

Explicitly, D(H) may be expressed in terms of the boundary values g(@) 
of H on the unit circumference c by the formula 


A(g) = — if dédg. (2) 


4 


Since, for given boundary values, a harmonic function gives the least 
value to Dirichlet’s integral, the principle (1) is equivalent to 


= af fie + dv = 


where the vector r(u, v) may represent any piece-wise continuously dif- 
ferentiable surface bounded by I 

If n = 2, so that I is a plane curve enclosing a region R, then, as was ob- 
served by the author, the Plateau problem becomes the one of mapping R 
on C conformally. Denote by x(u, v), y(u, v) the components of r; then 


Cc 


is a constant, equal to the area of R, for all vectors x(u, v), y(u, v) in the 
unit circle which map the circumference ¢ topologically on I’. Hence the 
principle (3) is equivalent to : 
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D(r) — A = min., (5) 


or 


E(x, y) -ff. [(x, — yp)? + (x, + y,)?] du dv = min. (6) 
Indeed, for a conformal map the Cauchy-Riemann equations 
x, — ¥, = 0, x, +y, = 0, (7) 


are satisfied, and E(x, y) = 0; while for a non-conformal map, E(x, y) > 0. 
The resemblance which the minimum principle (6) bears to the one 
considered by Riemann in his dissertation :* 


is only superficial. But recently the statement has been made, in refer- 
ence to (6), that: ‘‘This is exactly the famous variational problem which 
Riemann considered in his doctor’s thesis.’’ 

In view of the classic character of Riemann’s dissertation, and its influ- 
ence on the entire subsequent development of the theory of conformal 
mapping, it therefore seems appropriate to point out a number of sharp and 
essential differences between the principles (6) and (8). To do this is the 
purpose of the present note. 

2. Riemann was seeking a minimum principle to use as a basis for estab- 
lishing the existence of Green’s function G(P, A) of the region R, or the 
equivalent conformal map of R on the circular region C. 


The Green’s function is uniquely defined by the properties: (1) harmonic 
in R; (2) vanishing on T; (3) regular in the interior of R, except for a 
logarithmic singularity at the point A: 


G(P, A) = log p + h(P, A), (9) 


298 
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where p = distance AP, and h is a harmonic function of the point P regular 
throughout R including the point A. If G, denote the conjugate harmonic 


function to G, i.e., 
P 
PO 


er oe conformal transformation from R to C is given by the function 
e +34 

If [ is a general type of boundary curve, then the behavior of 0G/dn, 
therefore of G;, as we approach the boundary is uncertain. Hence the con- 
formal map furnished by Riemann’s method subsists only between the 
interiors of the regions involved, and the problem of what type of corre- 
spondence is established between the boundaries is left open. This problem 
was solved many years after the work of Riemann by Osgood and by Cara- 
théodory,® who proved that the conformal map of the interiors induces by 
continuity a one-one continuous correspondence between the boundaries. 

3. The first suggestion for determining G by means of a minimum condi- 
tion might be to use the Dirichlet principle as it stands: 


D(w) = ce (w,? + w,?) dxdy = min., (11) 


where we admit to comparison all functions w(x, y) of like nature to G(P, A): 
vanishing on the boundary I, and logarithmically singular at A. But since, 
as is easily verified, 


D(log p) = + © (12) 
over any neighborhood of A, we have the fundamental difficulty that 
D(w) = +o (13) 


for all functions w logarithmically singular at A. Hence the original 
problem (11) has no sense. 

We need a functional that for at least one value of its argument takes a 
finite value. Riemann secures this by means of an ingenious device, in 
the use of which (or similar devices) he has been followed by many subse- 
quent contributors to the theory.°® 

He introduces, instead of D(w) = min., the problem (8) 2 (w) = min., 
where 8 is a fixed function so designed that Q (w) is finite for at least a par- 
ticular form of w. Regardless of the parameter function B, this problem has 
the same variational equation as D(w) = min., namely, 


== (w, ia B,) + — (w, + B,) = Wr, + Oy = 0, (14) 
ox oy 


: 
= 
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or Laplace’s equation. The function 8 simply cancels out in the formation 
of the variational condition. 

Riemann’s idea was to arrange 8 so that Q(w) will be finite for at least 
one particular form of w. This he accomplished by choosing 6 to coincide 
in a small circle’ K about A as center with the conjugate harmonic func- 
tion to the one, log p, that characterizes the singularity at A, namely, 


B(x, y) = tan-"= (15) 


(we have taken A to be the origin). In R — K, 8 may be any smooth func- 
tion which (1) like tan~! y/x, has a period of 27 with respect toa cut Lin R 
that leads from A to the boundary I, and (2) which attaches smoothly to 
tan—! y/x at the circumference of K. The simplest choice is to define 8 by 
the formula (15) throughout the region R as provided with the cut L. That 
8 is multiform is no serious inconvenience, since only its partial derivatives 
are used, which are uniform: 


B, B, = (16) 


Riemann then selects as particular form of w a function a defined and 
smooth throughout R except at the point A, equal to log p in the circle K, 
and reducing to zero on I.’ Then, precisely because of the prearranged 
conjugacy between § and log p in the circle K, we have in that circle. 


Qa, — B, = 0, ay + B, = 0; (17) 


therefore 
Q(a) = 0 over K. (18) 


Since the partial derivatives of a and 8 are bounded and continuous? 
in the rest of the region, R — K, it follows that 


Q(a) = finite quantity. (19) 


Thus Riemann concludes that the problem (8) has sense. The minimiz- 
ing function w*, if existent,!° will then obey Laplace’s equation, and this, 
together with its vanishing on I and logarithmic singularity at A, charac- 
terizes w* as Green's function G(P, A). 

4. The preceding discussion may be summarized as follows. 

In the problem Q(w) = min. considered by Riemann, 8 1s entirely fixed, 
and the only function in competition is w, which varies subject to the fixed 
boundary values zero and logarithmic singularity at A. 

But in the problem E(x, y) = min., 

(1) x(u, v), y(u, v) are both variable in competition with other function- 
pairs; 
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(2) the boundary values of x, y are not fixed, but vary simultaneously 
subject to the restriction of always defining a topological correspondence 
between I and c. 

More important is the following distinction. 

(3) In Riemann’s problem, the argument function w(x, y) is defined in 
the generally complicated region R, and the sole and entire purpose of the 
problem is to obtain the Green’s function of R by identifying it with the 
minimizing function w”*. 

But in the problem E(x, y) = min., the argument functions are defined in 
the circular region C, where we have no need to find Green’s function, since 
that is known by elementary geometry: 


(20) 


Ay denoting the inverse point of A, and J the point (1, 0). 

Since E(x, y) differs only by a constant from the sum of Dirichlet integrals 
D(x) + D(y), the minimizing functions x(u, v), y(u, v) must be harmonic, 
and are therefore determined by Poisson’s integral as soon as their boundary 
values x(0), y(@) are known. Thus, the real point in connection with 
E(x, y) = min. is to find what topological correspondence shall be estab- 
lished by x(6), y(0) between I and c. 

The author has expressed this in an explicit way by the principle (2) 
A(g) = min., where g denotes the vector with the components x(6), y(@). 
The direct treatment of this minimum problem is safe from the Weierstrass 
criticism,!! since the argument range [g] is compact,'* whereas the argu- 
ment range [w] of Riemann’s problem is no/.* Accordingly, the existence 
of the minimizing vector x*(@), y*(0) is easily established, and the conformal 
transformation from C to R is then given immediately by the Poisson inte- 
grals for x*(@) and y*(@), or by the equivalent Cauchy formula 


i + iy") =e). 
Sc 


(4) The minimum problem of Riemann gives, when solved, only a cor- 
respondence between the interiors of Rand C. Supplementary considera- 
tions of quite different character were necessary to extend this correspon- 
dence by continuity to the boundaries. In fact, the topological nature of 
the boundary correspondence was only established by Osgood and by 
Carathéodory a number of years after the rigorous solution of the Riemann- 
Dirichlet minimum problem by Hilbert.'* 

But in the problem E(x, y) = min., the very nature of the argument 
functions, as prescribed, must automatically give a topological correspon- 


G(P, A) = log —— — log 

AvP AI 


302 MATHEMATICS: J. DOUGLAS Proc. N. A. S. 


dence between the boundaries simultaneously with the conformal map of 
the interiors. 

5. In conclusion, no basis remains for identifying the problem (6) 
E(x, y) = min. with the one considered by Riemann in his famous disserta- 
tion. 


1 J. Douglas, Trans. Amer. Math. Soc., 33, 263-321 (1931); Bull. Amer. Math. Soc., 
39, 227-251 (1933). 

2 This is regardless of whether the locus of the point (x, y) covers the region R simply 
or multiply, or extends outside that region. To assure the existence of the curvilinear 
integral in (4), we may suppose T to be rectifiable, but the fact of the constancy of % 
applies when I’ is any Jordan curve. 

3B. Riemann, Grundlagen fiir eine allgemeine Theorie der Functionen einer verdnder- 
lichen complexen Grésse, Gottingen, 1851. Gesammelte Werke, 1892. See especially 
§§16, 18, 21. 

4R. Courant, Ann. Math., 38, 696 (1937), where the transformation from (3) to (6) 
is given. See also Courant-Hilbert, Methoden der Mathematischen Physik, 2, 536-537 
(1937). 

5 W. F. Osgood and E. H. Taylor, Trans. Amer. Math. Soc., 14, 277-298 (1913); 
C. Carathéodory, Math. Annalen, 73, 305-320 (1913). 

6 Notably Weyl; see Die Idee der Riemannschen Fliche, 1913, p. 92. 

7 One whose radius is less than the distance from A to. 

8 In Riemann’s presentation, a + 18 is chosen equal to log (x + ty) in the circle K 
(called 6 by Riemann), which is equivalent to the separate conditions a = log p, B = 
tan—! y/x. R‘emann’s stipulation (loc cit., p. 41) that a + 78 be pure imaginary on the 
boundary of R is, of course, equivalent to: a = 0 on T. The region R is called T by 
Riemann. 

® We may even take a to be zero in a region adjacent to the boundary I’, as well as 
upon it. 

10 We have in mind the remark of Weierstrass, which emphasized the necessity of 
proving rigorously the attainment of the minimum. This was not accomplished until 
about 1900 by Hilbert, ‘‘Uber das Dirichletsche Prinzip,’’ Math. Annalen, 59, 161-186. 

11 See footnote 10. 

12 Every sequence of elements g contains a convergent sub-sequence. 

'3 For a fuller discussion of this point, see the second reference in footnote 1. 

14 See footnotes 5 and 10. 
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